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O : ABSTRACT 

We measured temperatures, gravities, and masses for a large sample of blue 
K> , horizontal branch stars in u Centauri, comparing the results with theoretical ex- 

pectations for canonical and He-enriched stars, and with previous measurements 
in three other clusters. The measured gravities of u Cen stars are systematically 
lower than canonical models, in agreement with expectations for He-enriched 
stars, and contrary to that observed in the comparison clusters. However, the 
derived masses are unrealistically too low as well. This cannot be explained by 
low gravities alone, nor by any of the other parameters entering in the calcu- 
lation. We find that the same stars are not brighter than their analogs in the 
other clusters, contrary to the expectations of the He-enrichment scenario. The 
interpretation of the results is thus not straightforward, but they reveal an intrin- 
sic, physical difference between HB stars in u Cen and in the three comparison 
clusters. 
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Subject headings: Stars: horizontal-branch — Globular clusters: individual (NGC 5139) 
— Stars: evolution — Stars: atmospheres — Stars: fundamental parameters 



INTRODUCTION 



The nature of the second parameter, aside from metallicity, which determines the mor- 
phology of the horizontal branch (HB) in globular clusters (GCs), is one of the most long- 
standing problems of modern astrophysics. In fact, a lower metallicity favors the formation 
of hotter and blu er HB stars, but clusters with the same metallicity can show very different 
HB morphology (ISandage fc Wildeylll967l), and an HB ex tended far toward the blue is ob- 
served even in some metal-rich GCs (e.g.. lRich et al.lll997l ). The helium abundance has been 



early pro posed, among others, as this second parameter ( ISweigartI Il997l : iD'Antona et al. 



2002 



see 



CatelanI l2009l for a review) because, during the He-burning phase, helium-rich 



stars are expected to be hotter than objects of canonical composition. This model has re- 
cently di;aw2_iMch_atteritioii, triggered by the disco v ery o f multiple stellar populations in 
GCs dPiottoet allboosi boOTh . In fact, IPiotto et al.l J2005I ) showed tha t a different met a l- 
licity is not the cause of the main-sequence split observed in u Centauri (JBedin et al.ll2004l ). 
and the only explanation is that the bluer sequence is greatly e nriched in helium, about 



50% more He-rich (Y=0.38) than in normal metal-poor GC stars (JNorrisll2004J : iPiotto et al. 



20051 ). In this scenario, the blue HB stars observed in many GCs could be the progeny 
of the He-enriched second stellar generation. Unfortunately, diffusion proces ses completely 
alter the surface chemical abundances of hot HB stars (e.g.. lBehr et al.lll999l ). preventing a 
direct demonstration of the connection between multiple populations and HB morphology. 
Nevertheless, an increased helium content can be indirectly deduced fr om other obser vable 
quantities, because He-enriched HB stars are predicted to be brighter (ISweigartI 119871 ) and 



to occupy different loci in the temperature-gravity plane ( Moehler et al.ll2003l ). 



In this Letter, we present the results of our investigation aimed to search for an in- 
direct indication of helium enrichment among blue HB stars in u Centauri, to test the He- 
enrichment scenario and its predicted effects on the HB morphology. This cluster is the ideal 
target for our purpose, because it hosts a very complex ste llar population, comprising three 
known MSs and six sub-giant branches (JBellini et al.l 120101 ). 



^Based on observations with the ESO Very Large Telescope at Paranal Observatory, Chile (proposal ID 
076.D-0810) 
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OBSERVATIONS AND DATA ANALYSIS 



We selected 116 target stars from the ground-based photometry of iBeUini et aLl ( 120091 ) . 
They span a wide portion of the cluster HB, from the blue edge of the RR Lyrae gap to 
the Blue Hook objects (Teff > 33 000 K) 5 mag fainter. In this Letter, we will focus on 
the comparison of uj Cen stars with other clusters and theoretical models. We therefore 
limit the analysis to Teff <33 000 K, because hotter Blue Hook stars are not included in 
the canonical models and are not present in the comparison clusters. A full analysis of the 
results, including the Blue Hook, will be presented in a forthcoming paper. 

The data were collected at Paranal Observatory in service mode between 2006 January 
and April, with F0RS2@UT1 in MXU mode. The selected 600B grism, coupled with O'.'S- 
wide slits, gave spectra of resolution R^ 1600 from the atmospheric cutoff to approximately 
5900 A. Three 45 minute spectra for faint stars and two 45 minute spectra for bright ones 
were acquired. Data were reduced with the FORS pipeline!, and the spectra were extracted 
under IRAFg, subtracting the nearby sky spectrurn with in the same slit. Finally, the spec- 
trum of the standard star LTT4816 (JHamuy et al.lll992l ). secured during observations, was 
used to flux-calibrate the object spectra, whose resulting signal-to-noise ratio was between 
50 and 150. Heli ocentric radial velocit ies (RV) were measured with the IRAF fxcor task, 
cross-correlating ( jTonry fc David 119791 ) the spectra with synthetic templates of adequate 
parameters as estimated from the stellar position in the color-magnitu de diagram (CMD) . 



Considering the internal velocity dispersion of the cluster (^13kms~^. ISoUima et al.ll2005[ ) 
and the errors of measurements (about 30 km s~^), the RV of all the observed stars is 
consistent with cluster membership. 

The atmospheric parameters of target stars were measured fitting th e observed Ba lmer 
and helium lines with stellar model atmospheres, computed with ATLAS9 ( lKuruczlll993l ). We 
used Lemke's versioru of the LINFOR program (developed originally by Holweger, Steffen, 
and Steenbock at Kiel University) to compute a grid of synthetic spectra. Stars showing 



iron lines in the 4450-4600 A region, indicating active atmospheric diffusion ( jMoehler et al. 



19991 ). or being hotter than 13 000 K (as deduced from the position in the CMD), were fitted 
with metal-rich models ([M/H]=-|-0.5) with varying surface helium abundance, to account 



2 www. eso.org/sci/data-processing/software/pipelines/index. html 

■^IRAF is distributed by the National Optical Astronomy Observatories, which is operated by 
the Association of Universities for Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 



For a description see http://a400.sternwarte.uni-erlangen.de/~ai26/linfit/linfor.html 
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Fig. 1. — Upper panel: distribution of uCen stars in the temperature-gravity plane. The 
Zero- Age and Terminal- Age HB (ZAHB and TAHB, respectively), for both canonical and 
He-enriched models are also indicated. Lower panel: comparison of stars in u Cen (full dots) 
and members of three other clusters (open circles). The vertical coordinate is the difference 
between the stellar gravity and the corresponding value of the canonical ZAHB at the same 
temperature. The plot is thus analogous to the Tcg-log (g) space of the upper panel, but 
the horizontal axis coincides with the canonical ZAHB. 
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for the effects of radiative levitation of lieavy elements (e.g., iMoeliler et al.ll2000l ). Tliis 
was done even for five stars hotter than 11500 K not satisfying these criteria, because the 
observed helium lines were clearly too weak compared to the fitted models when diffusion 
was not taken into account. The cooler stars were fitted with cluster- metallicity models 
([M/H] =— 1.5) with helium abundance fixed at the solar value, because their He lines are 
very weak and not observed at our resolution. In a few cases of doubt about the correct model 
to use, we adopted the set of model spectra that returned the lower y^ of the fit. The best fit 



to the obse rved spectra was ob tained by mean s of the routines develope d by [Bergeron et al. 



(I1992I ) and ISaffer et al.l (119941 ). as modified by iNapiwotzki et al.l (Il999l ). which employ a y^ 
test. The spectral lines used in the procedure included the Balmer series from H^ to H12, 
except for H, to avoid the blended Call H line, four He I lines (4026 A, 4388 A, 4471 A, and 
4921 A) for stars whose helium abundance was not kept fixed, and the Hell lines 4542 A 
and 4686 A when visible in the spectra of the hottest s tars. The routines es timate the errors 
on the derived parameters from the x^ of the fit (see iMoehler et al.l 119991 ) , but neglect all 
other sources of errors (e.g., defects in normalization, flat-field correction, sky subtraction). 
We therefore obtained a better estimate of the true errors multiplying the output values by 
3 (R. Napiwotzki, private communication). 

Stellar masses were calculated from the measured temperatures and gravities, through 
the relation: 



log 



M 



log 4 ■ log — — h log -— 



9q 



T^ 



^0 



where 



log 



U 



-0.4 -(V^ 



\m 



M)o - 3.1 ■ E[B -V)^BC- 4.74). 



(1) 



(2) 







We assumed T^=577 7 K, loggrT.= 4.44, (m-M)o=13.75±0.13 (Ivan de Yen et all 120061 ). and 
E(i? — V)=0.12±0.01 (JHarrislll996l . 2010 December Web version). The bolometric correction 
(BC) was derived from effective temperatures through the empirical calibration of iFlower 
(119961 ). Errors on masses were derived from propagation of errors. 



3. RESULTS 



Our results are plotted in the upper panel of Figure [11 where we show the position of 
the program stars in the temperature-gravity space, superimpos ed to the theor e tical zero- 
age and terminal-age HB (ZAHB and TAHB, respectively) from iMoehler et al.l (120031 ) . for 
canonical 
stars from 



Y=0.23) and He-enr iched (Y=0.33) models. In the same figure we include 78 



Moehler et al.l (120111 ). who measured the atmospheric parameters with the same 



procedure as in the present work, but based on medium-resolution FLAMES spectra, and 
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a different set of model spectra for stars above 20 000 K. Tlie two sets of data beliave very 
similarly in the Tcfj-log (g) plane. The comparison of the eight stars in common confirms 
the good agreement between the two works: the mean difference in gravity is null (< 0.01 
dex), while the difference in temperature is small (155 K), and becomes negligible (25 K) 
after the exclusion of two stars with very large errors (>1500 K). Our mass estimates are 
on average higher by 0.035 M (7), an offset accounted fo r by the fainter magnitudes (0.09 mag 
the mean difference) from the ICastellani et al.l (120071 ) catalog used by Moehler et al. Given 
the excellent agreement, we will merge the two datasets together. 

The surface gravity of stars cooler than ~18 000 K is systematically lower by about 0.2- 
0.3 dex with respect to canonical models, while they closely follow the trend of He-enhanced 
models at all temperatures. In the lower panel of Figure [T] we co mpare these results wit h 
similar measurement s obtained in three GC s, namely NGC6752 (JMoni Bidin et al.l 120071 ). 
M80 and NGC5986 (JMoni Bidin et al.ll2009l ). We adopt as vertical coordinate the difference 
between the measured log (g) and the value of the canonical ZAHB at the corresponding 
temperature. The comparison reveals that the HB stars in u Gen and in the other GGs behave 
very differently. We note that the stars of the comparison clusters do not completely agree 
with canonical models, whose expectation is to find the majority of the objects next to the 
ZAHB, and not to the TAHB as observed. Even so, u Gen stars clearly show lower gravities 
- at a given effective temperature - with respect to stars in other GGs. Observational 
errors tend to mask the general trend, but there is an offset of 0.15 dex at the cooler end, 
which decreases at higher temperatures and fades out around 18000 K, to reappear even 
larger (>0.2 dex) among hot stars at 25 000-28 000 K. However, we find a problem in the 
estimate of stellar masses, summarized in the upper panel of Figure El while the results in 



the comparison clusters roughly agree with expectations (see iMoni Bidin et al.l 120071 . 12009 



for a complete discussion), the mass es of uj Gen stars are const antly underestimated at all 
temperatures. Interestingly enough, iMoehler fc Sweigartl ( l2006l ) found very similar results 
for HB stars in N GG6388. This is a very peculiar cluster, as it has an extended, blue HB 
( iRich et al.lll997l ) despite its high metallicity ([Fe/ H]=— 0.6). This has been interpreted in 
terms of an extreme He-enrichment (up to Y=0.40. iGaloi fc D'Antonal 120071 ) , as in the case 
of its twin cluster, NGG6441. While iMoehler fc Sweigartl ( 120061 ) cast doubts on their results 
due to the large uncertainties caused by stellar crowding in this compact, bulge cluster, these 
results are very similar to what we find now in u Gen. 
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Fig. 2. — Upper panel: stellar masses, calculated through Equation ([T]). The line indicates 
the canonical model expectation. Lower panel: a bsolute magnitudes, est i mated from dis- 
tance moduli and reddening given in Section [2] and iMoni Bidin et al.l ( 120071 . 120091 ). Different 
symbols are used for u Cen and other cluster stars, as in the lower panel of Figure [H 



DISCUSSION 



In the temperature-gravity plane, the u Cen HB stars match the expectations of He- 
enriched models rather than canonical ones. However, the resulting underestimate of their 
mass prevents us from straightforwardly concluding that this is evidence of helium enrich- 
ment. In fact, the progenies of He-rich stars ii i the HB phase should ri ot be noticeably less 
massive than stars of canonical composition (JD'Antona fc Caloill2004j). an d the difference 



at any temperature is expected to be tiny (<0.03 M©, iMoehler et al. 



20031). Moreover, the 



derived masses are on average well below the value required to ignite helium in the core 
(-0.45 Mq). 

The easiest interpretation of our observations is the presence of a systematic error, bias- 
ing the results toward lower gravities and, as a consequence, lower masses. However, in thi s 
work we used the same instrument, software, and models as lMoni Bidin et al.l ( 120071 . |2009| ). 
finding a cl ear difference b e tween a; Cen and the other clusters, while our measurements well 
agree with iMoehler et al.l ( 120111 ). who also investigated cuCen but with a different instru- 
ment, higher resolution, different models for stars hotter than 20 000 K, and only a subset 
of Balmer lines. Therefore, even if we cannot completely exclude an observational bias with 
respect to theoretical expectations, the difference between u Cen and the three comparison 
clusters must be real: HB stars in u Cen are intrinsically different to their analogs in other 
GCs. The same conclusion can be drawn even if the offset is a product of the inadequacy 
of the employed models: in this case, they would be reproducing sufficiently well the atmo- 
spheric structure of the HB objects in the comparison clusters, but not in u Cen, hence a 
physical difference would be present. 

The observed trends cannot be due to a wrong (hotter) temperature scale. In fact, for 
each star we translated the measured temperature in a reddening estimate, compar ir ig the 
observed {B — V) color to the theoretical value obtained interpolating the iKuruczl (119931 ) 
grid, for the same metallicities as the model spectra used in the fits. The average value is 
E(i? — l^)=0.114, in perfect agreement with the literature, and with no significant trend 
along the HB. A temperature scale hotter by 10% (5%) would have caused an overestimate 
of reddening by 0.04 (0.02) mag at 10 000 K. In Figure [3] we compare the Yip line profiles 
of two stars with similar temperature in M 80 and u Cen: the line core and depth are very 
similar, indicating no noticeably difference in temperature, but the star in M80, whose 
measured gravity is 0.40 dex higher, shows wider wings. This comparison indicates that the 
peculiar stellar gravities refiect a real difference in the spectra of the target stars. We are 
aware that some stellar parameters unaccounted for in our study, such as stellar wind and 
rotation, can cause wider line wings mimicking a difference in gravity, and this degeneracy 
cannot be avoided at our resolution. The underestimate of gravities can also be caused by 
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the stratification of elemental abundances due to diffusion processes (JLeBlanc et al.ll2010l ). 
However, a systematic difference in one of these parameters, affecting uj Cen but not the 
other clusters, would be even more puzzling, and it is harder to postulate at this stage. It 
must also be noted that the mass underestimate does not completely follow the trend of the 
difference in gravity: while this seems to vary with temperature as noted before, the masses 
in a; Cen are constantly underestimated by ~0.15 M© for Teg >10 000 K. For example, at 
Tgff ~16 000-18 000 K there is a clear offset in masses, but not in gravities. This indicates 
that low gravities may play a role, but at least another effect should be at work, causing the 
observed mass underestimate. 



All the comparison clusters have similar metalli city ([Fe/H]?^ — 1.5 lHarrislll996[ ). while 
a; Cen shows a large spread up to [Fe/H]=— 0.6 (e.g. JSoUima et al.ll2005l ). Nevertheless, it is 
unlikely that the higher metallicity is the origin of the peculiar results in u Cen: the largest 
differences are found for stars hotter than ~11500K, whose surface abundances are altered 



by diffusion processes. iBehrl ( l2003l ) and iPace et al.l (120061 ) showed that, in the presence 
of diffusion, the surface abundance patterns are very similar in clusters of very different 
initial metallicity. The stars in all the clusters should therefore show the same behavior 
independently of their primordia l metal content, especiall y at 15 000-16 000 K where d if fusion 
reaches its maximum strength ( Moni Bidin et al.ll2009l ). Moreover, iMoehler et al.l (2000) 
found no peculiarity in the measured gravity and mass of two stars in 47Tuc ([Fe/H] =— 0.7), 
although using low-metallicity models. We repeated the measurements assuming different 
values of the model metallicity, to test how this parameter can affect the results. We found 
small differences in the stellar parameters, but the general behavior was unaltered: a higher 
model metallicity indeed returned slightly higher gravities, but higher temperatures too. 
As a consequence, the points were shifted almost parallel to the theoretical tracks in the 
Toff-log {g) plane, while the masses were increased by less than 0.05 Mq. 

The blanketing effect should be lower in metal-poor stars than in the solar-metallicity 
stars used to calibrate the adopted Teff-BC relation, and this could cause the underestimate of 
the BC and of the mass. T he adopted BC should be a good approximation for the stars hotter 
than the Grundahl jump (JGrundahl et al.lll999l ). because radiative levitation increases their 
surface abundances to super-solar values. As already noted, the effect should be independent 
of their primordial metallicity, thus the BC does not explain the offset of u Cen with respect 
to the other GCs. At cooler temperatures, the offset could be explained by the BC if u Cen 
stars were more metal-poor than in the other GCs, thus decreasing th eir BC by ~0.4 mag . 



Indeed, a; Cen hosts a metal-poor sub-population at [Fe/H]=— 2 (e.g.. iPancino et al.ll201ll ). 
but the BC varies by less than 0.15 mag for stars at 8 000-10 000 K in the w hole range between 
solar metallicity and [Fe/H] =— 2 (jCassisi et al.lll999l : lAlonso et al.lll999[ ). A wrong distance 



modulus or reddening could also cause wrong mass estimates, but the required correction 
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to increase the masses by 0.1 M© is huge (>0.4 mag in distance modulus, >0.13 mag in 
reddening). The recent hterature estimates agree on these quantities within 0.1 and 0.01 
mag, respectively, allowing only negligible variations of the mass estimates. Even an offset 
in the photo n aetric zero-point could c ause the observed offset , but comparing the V ma g of 
BeUini et all J2009[ l with the values of ICastellani et all feoOTJ l and iMomanv et all ^M) we 
found only a small difference (<0.1 mag) in the direction opposite to that required, with the 
magnitudes used here being brighter than in the other catalogs. In conclusion, none of the 
photometric parameters (stellar magnitude, BC, distance modulus, and reddening) entering 
in the calculation of the stellar mass through Equation ([2]) offers a viable explanation of our 
results. 

He-enriched HE stars are expected to be brighter than a nalogous objects of canonical 
composition (e.g.. lCaloi fc D'Antonall2007tlCatelan et al.ll2009[ ). and the increased luminosity 
should balance the lower gravities in Equation ([1]), returning a similar mass. All other 
parameters being the same. My should be brighter by 0.25-0.38 mag to compensate a decrease 
of 0.10-0.15 dex in gravity. On the contrary, in the lower panel of Figure |2]the perfect match 
between the absolute magnitudes of the HB of u Cen and the other clusters is clear. The too 
low values obtained for the stellar masses could therefore also be interpreted as due to the lack 
of increased flux, instead of too low gravity estimates. It could be argued that He-enriched 
stars could not necessarily be brighter in the V band, because the bolometric luminosity 
is the quantity entering in Equation ([1]). Thus, the detailed spectral energy distribution 
(SED) of He-enriched and canonical stars is required to properly deduce their luminosity 
from the V magnitude through Equation ([2]). However, great differences are not expected 
for stars hotter than 12 000 K, because the diffusion processes decrease the atmospheric He- 
abundance well below solar values in both cases. In fact, we find no difference in surface 
helium abundance between uCen and the other clusters, and \og{NHe/^H) < —1.5 for all 
the stars. With atmospheres of very similar chemical composition, th eir SED should not be 
very different, and even the known UV-enhanced flux of these stars ( jGrundahl et al.l 119991 ) 
should have the same effects irrespective of the primordial helium content. 



5. CONCLUSIONS 



Blue HB stars in uj Cen show lower gravities with respect to both canonical models and 
analogous stars in other GCs, but their stellar masses are underestimated, and their visual 
absolute magnitudes are very similar to those of the comparison clusters. Neither the low 
gravities nor the other parameters involved in the calculation can explain the too low masses. 
We can firmly conclude that these results reveal an intrinsic difference between the blue HB 
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stars in oo Cen and their analogs in other GCs, but its interpretation is not straightforward. 
The lower gravities follow the expectations for He-rich stars, but the magnitudes and masses 
do not. 
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Afines (CATA) BASAL PFB/06. G.P. acknowledges support by MIUR under the program 
PRIN2007 (prot. 20075TP5K9), and PRIN-INAF 2009. The authors are grateful to the 
anonymous referee for his/her fruitful report. 
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Fig. 3. — Comparison of the H^ lines of one of our program star (full line) and a HB star in 
M 80 at the same temperature. Derived temperatures and gravities are given. 



